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ABSTRACT 

One possible mechanism for the formation of molecular clouds is large scale colliding 
flows. In this paper, we examine whether clumpy, colliding, flows could be responsible 
for the clump mass functions that have been observed in several regions of embed- 
ded star formation, which have been shown to be described by a Salpeter type slope. 
The flows presented here, which comprise a population of initially identical clumps, 
are modelled using smoothed particle hydrodynamics (SPH) and calculations are per- 
formed with and without the inclusion of self-gravity. When the shock region is at its 
densest, we find that the clump mass spectrum is always well modelled by a Salpeter 
type slope. This is true regardless of whether the self-gravity is included in the simu- 
lations or not, and for our choice of filling factors for the clumpy flows (10, 20 & 40%), 
and Mach number (5, 10 & 20). In the non-self-gravitating simulations, this slope is 
retained at lower Mach numbers (M.= 5 & 10) as the simulations progress past the 
densest phase. In the simulations which include self-gravity, we find that low Mach 
number runs yield a flatter mass function after the densest phase. This is simply a 
result of increased coagulation due to gravitational collapse of the flows. In the high 
Mach number runs (A4— 20) the Salpeter slope is always lost. The self gravitating cal- 
culations also show that the sub-group of gravitationally bound clumps in which star 
formation occurs, always contain the most massive clumps in the population. Typi- 
cally these clumps have a mass of order of the Jeans mass of the initial clumps. The 
mass function of these bound star forming clumps is not at all similar to the Salpeter 
type mass function observed for stars in the field. We conclude that the clump mass 
function may not only have nothing to do with gravity, but also nothing to do with 
the star formation process and the resulting mass distribution of stars. This raises 
doubt over the claims that the clump mass function is the origin of the stellar IMF, 
for regions such as p Oph, Serpens and the Orion B cloud. 

Key words: stars: formation - ISM: clouds - ISM: kinematics and dynamics - ISM: 
structure 



1 INTRODUCTION 

The mass spectrum of stars, clumps or clouds, is normally 
described by some sort of power law, which has the form 
dN oc m 1 dm, where N is the number of objects, and m 
their mass. Recent sub-millimetre observations of clumps 
and cores in the star forming regions of p Ophiuchus, Ser- 
pens, and Orion B have shown that the mass distribution of 
these gas density enhan cements is similar to the mass dis- 
tribution of young stars (| Motte e t al.lll998l : iTesti fc Sargenj 
ll99Sl:IJohnstone et al.l200ll200ll) . More r ecent observations 
using the extinction mapping technique jAlves et alJl200lf) 
have confirmed this res ult in a number of star forming re- 
gions iLada et al ] l200dh These observations have been pre- 
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sented as evidence that the fragmentation process controls 
the mass of stars. In this picture, a bound region of smooth 
gas breaks up into a series of fragments (the clumps or cores, 
depending on your terminology) as the gravitational insta- 
bility sets in. These fragments then collapse to form individ- 
ual stars. 

There is however the suggestion that most clumps in 
this region appear to be stable against gravi tational collapse, 
at least in p Oph and the Orion B cloud iJohnstone et all 
l200d 1200 lh . We therefore do not know if the clumps that 
make up the mass distribution are bound objects, or if they 
are just transient features that will never collapse to form 
stars. It is possible that most of the density structures ob- 
served in these regions have little to do with the star forma- 
tion process. 

Such transient density features are commonly found in 
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simulations of turbulent gas, leading to the suggestion that 
the clump mass spectra, and thus potentially the stellar 
IMF, is con trolled by shock density e nhancements in tur- 
bulent gas iPadoan fc Nordlundl 12002'). Numerical studies 
of driven and freely decaying turbulence show that a clump 
mass spectrum similar to the the IMF can be formed via the 
shocks arising from the turbule nt flows (e.g. lKlessenll200 it 
iBallesteros-Paredes et alll2005l) . It is generally found that 
most of these clumps are unbound. IClark fc Bonnelll ll2005h 
have also argued that the clumpy structure, and thus the 
turbulent motions themselves, are essentially irrelevant to 
the star formation process, with the mean Jeans mass in the 
cloud controlling the star formation. 

It is still unclear however whether the 'turbulence' mod- 
elled in the above star formation simulations is an accu- 
rate description of the velocities inside molecular clouds. For 
example, molecular clouds, and the giant molecular clouds 
(GMCs) where most star formation occurs, may not be co- 
herent structures, but rather an ensemble of smaller clouds 
which have been acc umulated by some external mechanism 
iBonnell et al.ll20 05l) . At the point of a cloud's 'formation', 
the smaller clouds out of which it comprises will probably 
have some random motions. In this model, the internal struc- 
ture and velocity of a parent molecular cloud is then dictated 
by the state of the constituent clouds. Regions of high den- 
sity, such as the regions where stars can form, will therefore 
be the result of clump-clump collisions, or of larger collisions 
between coherent flows of clumpy material. 

The collisions between the internal structures of these 
molecular clouds (the original smaller clouds) hav e been 
sugge sted as a possible trigger for star formation iSmitbl 
ll980Tl . There are extensive studies of cloud-cloud collisions 
in the literature, involving numerical calculations of colli- 
sions between pairs of clouds. These have looked at dif- 
ferent Mach numbers, collisions between clouds of differ- 
ent sizes, gas cooling during the sh ocks, and the appear- 
ance of young emb edded protostars (IChapman et al.lll992t 
iBhattal et al]IT998h . The results of these simulations sug- 
gest that is should be possible to form a wide range of binary 
separations and component masses via random collisions be- 
tween clumps of varying size /mas s etc. Only one study to 
date however. iGittins et al.l (I2003I) . has modelled an ensem- 
ble of such 'cloudlets' using a fluid code. The purpose of 
that study was to examine whether the coagulation process 
can lead to a stellar type clump IMF, assuming a system 
of initially identical cloudlets, with random motions. They 
found that very few mergers occurred in their simulations 
(due to the assumed small filling factor), resulting in a fi- 
nal mass spectrum of cloudlets that has a slope somewhat 
steeper than that of the Salpeter IMF. 

Cloud-cloud collisions have also been i nvestigated ana- 
lytically, via the coagulation equation (e.g. [Field fc Saslaw 
1 19651: iField fc Hutchins|ll96g| lKwanlll979l: IScoville fc Hersh 
Il979l : ISilk fc Takahashilll979T) . Assuming various combina- 
tions of cloud cross sections and velocities, these stud- 
ies concluded that the steady state mass distribu t ion o f 
the clouds would have 7 ~ 1.5. iPenston et al.l <ll969T) 
and lHandburv et alj Jl977t) . using Monte-Carlo simulations, 
found 7 in the range 1.5 to 1.9 respectively. However all these 
results assume that contact between two clouds will result 
in a m erger. Using 3D hydrodynamic simulations. iHausmanl 
dl98lT) demonstrated that not all collisions result in merg- 



ers, with high Mach n umbers causing the clouds to be shred- 
ded. [HausmiMi| lll982l) used these numerical collision calcula- 
tions to include a prescription of the shredding/merging pro- 
cess into Monte-Carlo simulations. Again, the results showed 
that the mass spe ctrum should have 7 in the range 1.5 to 2. 
lElmeereenl (Il989t) included the effects of self-gravity into the 
coagulation equation, and how it can increase the chances 
of cloud-cloud collisions. His results showed that the inclu- 
sion of 'gravitational focusing' of the collisions can result 
in a steady-state cloud mass spectrum with 7 = 2, but the 
calculations didn't include a prescription for possible shred- 
ding. In fact, by parameterising the various properties of the 
initial cloudlet popul ation, one can get a wide range of final 
mass function slopes JSilk fc Takahashilll979h . 

The results of both the analytic and the numerical 
studies have shown that a power law type mass spec- 
trum can be achieved from clump-clump (or cloud-cloud) 
collisions between objects moving with random motions. 
It has been suggested however t hat star forming regions 
are formed from large sc ale flows dVazouez-Semadeni et aJl 
Il995t IPringle et aljEoOlll . If this is true, it is unlikely that 
these flows would be of a uniform density, but instead be 
clumpy. In this scenario, the chances of collisions are greatly 
increased. The resulting mass spectrum is also unlikely to be 
steady state, as is found for the random motion studies, since 
the flows have finite extent. 

In this paper, we wish to examine whether the shock 
compressed regions created by supersonic, and clumpy, flows 
can produce a clump/core mass function that resembles the 
stellar IMF. We wish also to examine the importance of 
gravity in building the clump/core population. We do not 
attempt here to fully explore the potentially vast param- 
eter space associated with clumpy flows but instead con- 
sider the simplest case of initially identical, uniform density, 
clumps. We thus only try to test the concept here, rather 
than build up a detailed picture of cloud formation. Both 
self-gravitating and non self-gravitating simulations are pre- 
sented. In section|2]we present the details of the simulations, 
including a discussion of the SPH code and the method of 
clump-finding that we have used. We go on to discuss the 
results of the non self-gravitating calculations in Section |H] 
and the results of the simulations with self-gravity in Sec- 
tion 0] There is a discussion of these results in relation to 
the p Oph star forming region in Section |5] and we conclude 
with a brief summary of our main findings in Section [S] 

2 THE SIMULATIONS 
2.1 The fluid code 

The fluid is modelled using the Lagrangian part icle method 
of smoothed particle hydrodynamics, or SP H llLucvl Il977l : 
iGingold fc Monaghanlll977HMonaghanlll992l) . The smooth- 
ing lengths are variable in both time and space, with the 
constraint that there must be a roughly constant num- 
ber of neighbours for each particle, which is chosen to 
be roughly 50 (with a fluctuation from 30 to 70 neigh- 
bo urs). We use the standard ar tificial viscosity suggested 
by iMonaghan fc Gingoldl lll983T) with an alpha and beta 
of 1 and 2 respectively. Gravitational forces are calculated 
using quadrupole moments obtained via a tree structure 
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Figure 1. We plot here column density images taken from the Mach 10 simulations. The images are taken at a time t = 0, t = 0.5t C r, 
and t = t cr , where t cr is the crossing time for the flows, defined as t cr = L fiow/ v flow F° r these simulations, Lfi ow is just half the 
length of the region shown. The top, middle and bottom rows contain the images from the 10%, 20% and 40% filling factor simulations 
respectively. The colour table for these images is stretched over a range from 0.0001 g cm -2 to 10 g cm 



l|Benz et al.lll99(t l. which is also used to construct particle 
neighbour lists. 

The code in c ludes the modification by 
iBate. Bonnell fc Price! il995T) such that dense, bound, 
regions of gas are replaced by point masses, or 'sink 
particles'. The sinks interact only via gravity but are able 
to accrete gas particles within their accretion radius. These 
sinks allow the code to model the dynamical evolution of 
accreting protostars, without integration times becoming 
prohibitively small. Only three of the simulations presented 
here include gravitational forces and sink particles. The de- 
tail of the sink particle formation in these three simulations 
in given below. 



2.2 Initial conditions 

In this paper, we present a series of simulations following the 
evolution of colliding flows, in which the flows consist of a 
population of clumps. The clumps are created by populating 
the computational volume with a settled template clump of 
2000 particles. The clump radii are then set to provide a de- 
sired filling factor. Masses of the clumps are then determined 
by assuming each template clump has a Jeans number of J 
= 0.5 at a temperature of 10K. The Mach number of flows is 
also changed in the simulations, from Mach 5 to Mach 20. A 
summary of the initial conditions for each of the simulations 
can be found in Table [T] 



4 Clark & Bonnell 




Figure 2. Shown in this figure are the clump mass spectra for all the calculations without self-gravity. The spectra are shown at two 
instances during the flow evolution. On the left, the spectra are taken at a time t = 0.5t cr , corresponding to the point where the whole 
cloud is at its most compact (that is, the full length of the flows has entered the shock). On the right, we show the spectra at time 
t = t CT , which is the point where the flows have completely emerged from the shock. In each of the plots we represent the three filling 
factors. The solid, dotted and dashed distributions showing the spectra from the 10, 20 and 40% filling factor flows, respectively. The 
long dashed line shows a Salpeter slope of 7 = 2.35 (so —1.35 on our log-log plot here). 
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Table 1. Details of the nine simulations with no self-gravity 
arc presented below. The volume filling factor of the flow, /, if 
defined as the percentage of the flow volume that is occupied by 



M 



flow 



is the total mass in the flow, and M,, 



is the 



mass of the initial template clumps that are used to populate 
the flow's volume, and r c i ou di e t I s the radius of these clumps. 
The calculations which include self-gravity have identical initial 
conditions to those of simulations 7, 8 and 9. 



Sim No. 


/ 


M 




^cloudlet 


^cloudlet 




(percent) 




(M ) 


(M ) 


(pc) 


1 


10 


5 


1104 


4.4 


0.37 


2 


10 


10 


1104 


4.4 


0.37 


3 


10 


20 


1104 


1.1 


0.37 


1 


20 


5 


1390 


5.6 


0.46 


5 


20 


10 


1390 


5.6 


0.46 


6 


20 


20 


1390 


5.6 


0.46 


7 


10 


5 


1752 


7.0 


0.58 


8 


10 


10 


1752 


7.0 


0.58 


9 


10 


20 


1752 


7.0 


0.58 



We now go discuss this process of constructing the flows 
in more detail. The clumpy flows are formed by populating 
the computational volume with copies of a template clump. 
To construct the original template clump, we place 2000 
SPH particles randomly in a spherical volume. The particles 
are given no initial velocities. We then evolve this clump in 
isolation for ten dynamical times 1 with a spherical, constant 
volume, boundary condition. This allows the SPH particles 
to settle into a uniform density sphere and gives time for 
particle velocities, which result from the random initial po- 
sitions, to be damped by the artificial viscosity. When the 
template clump is used to create the clumpy flow, each copy 
is given a random orientation before being placed in the 
computational volume. 

Two properties of the flow are altered in this paper: the 
Mach number and the filling factor. The number of initial 
clumps is held constant at 250, and the volume of the flow is 
kept the same, with the computational volume being a cube 
of side 8 pc. To alter the filling factor, the volume of the 
clumps is then changed such that the clumps occupy the de- 
sired fraction of the computational volume. We also impose 
the constraint that the clumps are unbound by their inter- 
nal energy, with a Jeans number J = \E gra v\/ Etherm = 0.5. 
They clumps are also given a temperature of 10K, which 
remains constant throughout the simulation, since we as- 
sume an isothermal equation of state. The combination of 
filling factor and Jeans number then sets the required initial 
mass, and thus density of the clumps. Without some kind 
of external pressure, the clumps, being unbound, would just 
expand and eventually merge. To prevent this occurring, we 
set an external pressure term which is equal to the thermal 
pressure of the clump gas. 

The flow velocities are set up along the x axis, such that 
particles with x < have a velocity v x = +Vfi ow — +Mc s 
and particles with x > have a velocity v x = —Vfi ow — 
—Mcs- The y and z components of the velocity are set to 



1 The dynamical time is just the sound crossing time for the gas 
sphere. 



0. As mentioned above, the cubes have a side of 8 pc in 
our simulations, which goes from —4 to 4 pc in each of the 
x, y and z directions. In this paper we will refer to the 
crossing time, which we will denote by t CI . The crossing time 
is given by t cr — Lfi ow /vfi ow . Since our flows only take up 
half the length of the computational volume, the crossing 
time is t cr = 4 pc /(Mc s ). We define the crossing time in 
this manner so that the shock is finished at t = t CI . 

In this paper, three different flow velocities are pre- 
sented. These are Mach 5, 10 and 20, which with a sound 
speed of 0.2 km s - , gives flow speeds of 1, 2 and 4 km s - , 
respectively. The corresponding crossing times for the Mach 
numbers 5, 10 and 20 are 4, 2 and 1 Myrs, respectively. For 
each of these Mach numbers, we also simulate 3 different 
filling factors, with the clumps occupying 10, 20 and 40% 
of the flows. The filling factors will be denoted in this pa- 
per by the symbol /. The template clump which is used 
to populate the flows has a mass of 4.4 M©, 5.6 Mq and 
7 M for the 10%, 20% and 40% filling factor simulations, 
respectively. There are thus 9 different initial conditions for 
the simulations. In these simulations we do not include grav- 
itational forces, since our first aim is examine the ability of 
gas to form clump-mass spectra purely through collisions. 
However for the simulations discussed in sectional the self- 
gravity is switched on, allowing us to examine the effect of 
the gravity on the clump mass spectra, and the formation 
of stars. 

To follow the star formation in simulations which in- 
clude self-gravity, we create 'sink particles' when the den- 
sity reaches a certain limit. In SPH calculations involving 
self-gravity, there is a limiting mass resolution below which 
the S PH method will not r e liably form self-gr avitating ob- 
jects jBate fc BurkerlJll997l) . lBate et all (12003ft showed that 
this occurs when the objects are represented by less than 
roughly 75 SPH particles. The mass resolution of any SPH 
simulation can thus be given by, m reB ~ 100(Af s ; m /7Vj, ar t). 
In our calculations, we create the sink particles at the point 
where m res becomes self-gravitating, and thus has a Jeans 
mass. The Jeans mass is given by, 



mj 



-1/2 



3/2 



(1) 



5 /try 

2 G^J 

and one can rearrange for the density to obtain the den- 
sity threshold above which sink particles can be created. 
The mass resolution for our simulations which include self- 
gravity (with / = 40%) is ~0.3 Mq , which gives a sink par- 
ticle creation density of 2.47 x 10 -18 g cm -3 . The radius as- 
sociated with this mass resolution and density is 3.87 x 10 16 
cm (or 0.012 pc) and represents the accretion radius for the 
sink particles, inside which they can accrete gas particles. To 
ensure that our codes run quickly, we smooth the gravita- 
tional forces between sink particles to the accretion radius. 
This prevents binaries with separations smaller than 0.012 
pc from forming in our simulations, since they are not of 
interest to this paper. 

2.3 Finding the clumps 

The main aim of this paper is to examine the formation 
of clump populations in colliding flows. This study uses 
a clump finding algorithm based on the method devel- 
oped by IKlessen fc Burkertl i2000l) . This is in turn based 
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on th e method developed by I Williams, de Geus fc Blit3 
ll994|). We point the inte rested reader to the Appendix 1 in 
iKlessen fc Burkertl l)200(t) for a description of the method. 
This algorithm makes use of the SPH formalism and has 
been shown to very easily identify density structure in the 
gas. These structures are analysed in this paper in order to 
assess their role in the star fo rmation process . 

However unlike bo th IWilliams et alj (I1994T) and 
IKlessen fc Burkertl i2000l) . we do not use a density thresh- 
old to define what is termed as a 'clump', i.e. there is no 
condition on how much higher the peak density has to be in 
relation to the minimum density of the structure. While this 
means that, technically, very low contrast density features 
are free to be defined as clumps, in practise this never oc- 
curs in our simulations. The clumps, formed from supersonic 
flows, have a density peaks that are an order of magnitude 
higher than the mean gas density. 



3 THE CLUMP MASS SPECTRA 

In Figure Q we show column density images for the simula- 
tions in which the flow has a Mach number of 10. Each of 
the three filling factors that is investigating in this study are 
presented, and images are shown at three times during the 
simulations: first the initial flow setup, then at t = 0.5t cr (half 
crossing time, when the whole ensemble of gas is at its most 
compact), and at t =t CI , when the flows are roughly back 
to their original size. One can easily see that at the half 
crossing time, the cloud is compact, and that it is heav- 
ily structured with dense filaments that have formed as the 
colliding clumps produce shocks. After a crossing time, the 
cloud has less shock structure, and the density features are 
more rounded, since the shock filaments have had time to 
re-expand. The region now has a mottled and wispy ap- 
pearance, similar to what one observes in the gas structure 
around embedded clusters. 

The purpose of this study is to see whether the process 
of creating a dense region from a clumpy flow, can create a 
population of clumps in the region that resemble the mass 
spectrum for stars. The clump mass spectra are shown in 
Figure |21 This figure shows the mass spectra for each of the 
simulations, taken at t = 0.5£ cr and t =t CI . Note that the 
clumps shown in the figure, as found by the clump-finding 
algorithm (see Section [2.311 . are 3-dimensional. 

It is clear from the plots in Figure|2]that the clump mass 
spectra evolve as the flows enter the shock. The distributions 
at t = 0.5t cr and t =t CI look distinctly different. The amount 
of evolution in the spectra depends on the volume filling 
factor /, with the lower / simulations changing much less 
than those with higher values of /. This result is not that 
surprising. In the calculations with the high filling factors, 
there is a much greater chance of two clumps colliding and 
potentially shredding. 

The Mach number of the flows can also have an effect 
on the mass distribution of clumps. In the Mach 5 simula- 
tions, there are more high mass objects at t =i CI .than there 
are at t =t CI . The distributions are also better fitted by a 
power law. For the simulations with Mach 10 & 20 flows, 
the opposite is true. The more massive clumps appear at 
t —0.5t CI . The distributions is also better fitted by a power 



law at t = 0.5t cr , and by comparison, the clumps are less 
massive at t =t cr , and not easily fitted by a power law. 

This dependence on the Mach number for the clump 
mass spectra suggests that two different processes are at play 
here, coagulation and shredding, and that the efficiency of 
these processes is linked to the strength of the shocks. One 
thing that is certain is that the shredding is always much 
more efficient here than the coagulation. The initial clump 
that is used to populate the flows has a mass of roughly 
4 to 7Mq , depending on the filling factor of the flow in 
question, and there are only 250 clumps present initially, so 
the shredding process clearly dominates. 

At lower Mach numbers however (represented here by 
our Mach 5 simulations), coagulation does seem to be 
present. Low mass clumps gradually grow through merg- 
ing events with more massive clumps. The low mass clumps 
must therefore survive for long enough (that is, no re- 
expansion) to be able to keep merging with other clumps. 
For the clumps to get progressively more massive as the 
flows interact, implies that either the time scale for their re- 
expansion is less than t cr , or that the clumps are in equilib- 
rium with the pressure that we included in the simulations. 

At higher Mach numbers, there appears to be no such 
steady coagulation process. In the Mach 10 flows, the clump 
mass spectra don't appear to undergo much evolution be- 
tween times 0.5t C rand t CT . This suggests that the balance 
between shredding and coagulation, whatever it may be, is 
fairly constant throughout the simulations. It is interesting 
to note here that only the 40% filling factor flows manage 
to produce clumps that are more massive than those at the 
start of the simulation. When one moves to the Mach 20 
simulations, the situation is slightly different. At t = 0.5t cr , 
the mass spectra look quite similar to those in the Mach 
10 simulations, however the Mach 20 flows manage to form 
some more massive objects. At t =£ cr however, much of then 
higher mass structures have been lost. Either, these struc- 
tures were not in equilibrium with their surroundings, and 
re-expanded to become part of the surrounding objects, or 
they have been broken up by destructive encounters with 
other clumps. In either case, shredding and coagulation are 
not balanced in the Mach 20 simulations, and the higher 
mass features in the clump mass spectrum at t = 0.5t cr are 
only transient. 

As mentioned above, the purpose here is to determine 
whether the formation mechanism for forming regions of star 
formation is potentially the origin of the clump mass distri- 
butions within. We can see from from Figure |5J that many 
of the clump mass spectra can be fitted by a Salpeter type 
power law. It is therefore possible for colliding, clumpy, flows 
to produce a clump/core population with a mass spectrum 
consistent with that for the field stars. Clearly, from the 
spectra shown in Figure not all the simulations display 
this property at all times, but most do exhibit a Salpeter 
type slope at some point. 

It is important to note here that none of the resulting 
clumps will be the sites of star formation in these simula- 
tions, since the model contains no self-gravity, and gravity 
has played no part in producing the mass spectrum. All the 
structures represented in Figures Q and |2] are purely due to 
the coagulation and shredding processes that occur as the 
flows interact. 

These result suggest that the observed clump mass spec- 
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tra in a number of star forming regions of may just be a 
result of how the region was formed, and not directly part 
of the collapse and fragmentation process. 

Although many of the more massive objects are lost 
in the higher mach number flows, and the mass spectrum 
looks less like that for field stars by the time t =t CI , this is 
not a problem for the model. Star formation only occurs in 
the regions of highest density (> 10 5 cm -3 ) within molec- 
ular clouds, and, consequently, it is in these regions where 
the clump/core observations are made. In terms of our sim- 
ulations, the density is highest, and the cores much better 
defined, around the time t ~ 0.5t CI , and it is at this point in 
the evolution where the clump mass spectra are most similar 
to the stellar IMF. 



4 THE INFLUENCE OF GRAVITY 

In the previous section we showed that clumpy colliding 
flows can produce a region with a clump mass function that, 
in the absence of self-gravity, resembles the stellar IMF. We 
find that shredding dominates over coagulation in all these 
simulations, however in our lower Mach number runs (M = 
5) we see that some steady coagulation does occur and cre- 
ates a number of higher mass objects. In the higher Mach 
number runs {M= 10 & 20), higher mass clumps are also 
present at the point in the simulations where the cloud is at 
its most dense (t — 0.5t CI ). These are only transient features, 
however, and most are gone by t =t CI . 

In this section we examine the effects of self-gravity in 
some of these runs. We have only selected 3 simulations in 
this section, to be re-run with the self-gravity switched on, 
since the computational expense of this type of calculation 
is much greater. We have chosen to re-run the 40% filling 
factor simulations, with Mach numbers 5, 10 and 20. We 
have chosen these simulations since at Mach 10 and Mach 
20, their non-self-gravitating versions give mass spectra that 
are remarkably similar in slope, at t — 0.5t CI , to the stellar 
IMF. 

The three 40% filling factor simulations were re-run 
with the self-gravity switched on. Just like the original ver- 
sions, they were allowed to run for one crossing time. Figure 
13 shows column density images from the end of each of the 
three simulations. The Mach 5 and Mach 10 flows both re- 
sulted in star formation during the crossing time but no star 
formation occurred in the Mach 20 flow. In the simulation 
with the Mach 5 flow, the star formation occurred at t = 0.45 
t CI ( or 1.8 xlO 6 years), and in the Mach 10 simulation at 
t = 0.63 t CI (or 1.2 x 10 6 years). 

Figure [I] shows the clump mass spectra for these new 
simulations. At t = 0.5t CI the spectra look very similar to 
previous simulations where self-gravity was excluded. The 
Mach 10 and 20 flow simulations thus have the same Salpeter 
type slope, and now so too does the Mach 5 flow simulation. 
At t —t CT the Mach 5 and 10 simulations have an overall flat- 
ter distribution than before, having more high mass clumps, 
and fewer low mass clumps. The Mach 20 simulation has 
essentially the same distribution. 

The inclusion of self-gravity to the calculations therefore 
did not manage to allow the Mach 20 simulation to preserve 
the more massive clumps that are formed during the flow. 
However, its inclusion also did not significantly alter our 



previous result, that during the densest phase of the flows, 
that is time t = 0.5t CT , the mass spectra of the region has a 
Salpeter type slope. 

The fact that the clump mass spectra for the Mach 
5 and 10 flows is flatter at t =t cr than in the non-self- 
gravitating case is simply a result of the flows collapsing 
under self gravity. This collapse can be seen clearly in the 
column density images in figure |3] Instead of having a 1-D 
compression of the gas along the x axis, as occurred in the 
non-self-gravitating cases, the simulations with self-gravity 
experience a a 3-D compression. This results in more coag- 
ulation of smaller clumps into larger structures. The Mach 
20 simulation is less affected by this gravitational collapse, 
since the flow crossing time is much shorter than the gravi- 
tational free-fall time of the flows: it simply has less time to 
feel the effects of the self gravity. 

It is interesting that the Mach 5 clump distribution is 
also more Salpeter like with the inclusion of gravity. If we 
look back at the 10% filling factor runs in figure |2J we see 
that they are always steeper than those for the higher filling 
factor runs. This is probably a result of the higher initial 
clump density in the low filling factor runs (see section 
for details of how the runs are set up). At these higher den- 
sities, the shredding is likely to be more efficient, due to the 
increased the ram pressure. Similarly, in the 40% filling fac- 
tor run, with the self gravity included, the collapse of the 
flow under the gravity will create denser clumps, perhaps 
aiding the initial shredding, before the flow collapse can re- 
sult in coagulation. 

Figure 2] also plots the distribution of clumps for which 
the gravitational energy is greater than the thermal energy, 
which we call thermally bound, and the clump distribution 
for which the gravitational energy is greater than both the 
supporting kinetic and thermal energies. We will refer to 
these latter clumps as totally bound. We must remember here 
that none of the structures initially have any internal kinetic 
energy. The kinetic energy they now possess has come from 
the shock formation process. The thermally bound popula- 
tion generally spans a wider range in clump mass than the 
totally bound clumps. It is only the more massive clumps in 
the distribution that manage to have little enough kinetic 
energy to become totally bound. The star formation that 
occurs in these simulations is therefore confined to the most 
massive of the clumps. Note that these star forming clumps 
do not constitute a representative subset of the main clump 
distribution. In fact, a similar result was seen in simulations 
of bo th driven (iKlessenl l200lf) and freely decaying turbu- 
lence JClark fc BonneJl2005ft ~ in that in both cases the star 
formation was confined to largest clumps. 

The fact that the bound clumps form at or around 
20 M© shows that the shock clump-formation process does 
little to alter the energy state of the original clump gas. 
The clumps originally had a Jeans number of 0.5 Mq, 
which with an original clump mass of 6 Mq, gives the 
original Jeans mass in the clump gas to be 17 Mq. Thus 
only clumps that have a mass similar to the original Jeans 
ma ss become bound. Th is is consistent with the findings 
of JClark fc Bonnellll2005ft . who showed the same result for 
freely decaying turbulence. 

As already mentioned, the Mach 20 flow simulation 
failed to produce any stars in a time period of one cross- 
ing time. Interestingly, it did have one bound clump at 
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Figure 3. Shown here the column density images of the 40% filling factor simulations, this time re-run with the inclusion of self-gravity. 
The images are taken at the crossing time for the flow, the point at which the simulations were terminated. The maximum and minimum 
column densities plotted are 10 and 0.0001 g cm -2 , respectively. 



t — 0.5t CI , which (again) was among the most massive of 
the clumps in the simulation, but this clump was lost as 
the flow evolved. At the point at which the simulation was 
stopped, the clump population does have quite a few lower 
mass members that are thermally bound. It is possible that 
if the flow was left for longer that this simulation might 
eventually produce some stars. 



5 DISCUSSION 

We have proposed in this paper that the clump mass spectra 
that have been observed in p Oph, Serpens and the Orion B 
cloud may have arisen naturally if the regions were formed 
from large scale shocks, p Oph also display s some evidence 
for a shocked formation. iMotte et alJ Jl998f) discuss the pos- 
sibility that p Oph has been triggered b y a large shock from 
the nearby OB association of Sco OB2 dde Geuslll992l) . 

The morphology of p Oph alone is suggestive of a shock 
formation process. The region has a filamentary structure 
in the form of ridge of hig h column density gas, first shown 
bv lWilking fc Ladal II1983) from 18 CO observations. Higher 
density molecular tracers also revealed that this ridge is bro- 
ken up into a series of dense cores (labelled/ r eferred to as 
Oph-A to Oph-F.), lLoren fc Woottenl lll986l) : lLoren et all 
il990l) . These cores are also clearly visible in the 1.3 m m 
continuum observations presented bv lMotte et alJ Jl998l) . 

It also appears that the region, and the star formation 
within, occurred fairly rapidly in p Oph. The young stellar 
objects, buried in the cores found along the ridge, have a rel- 
atively strong 1.3 mm continuum emission, suggesting that 
they are young Class I protostars, with ages of 10 5 years. 
This is an order of magnitude less than the crossing time 
along the ridge, suggesting that the star formation occurred 
in a fairly rapid burst. The star formation that occurs in 
the simulations we present in this paper also occurs fairly 
rapidly, on a timescale less that the flow crossing time. Both 
the observations of p Oph and the simulation s presented 
here a re therefore consistent with suggestion by lElmegreenl 
teOOCf) . that star formation occurs on a timescale compara- 
ble to the crossing time in the region. 

Although the triggering for the p Oph region has so 



far only been attributed to the nearby Sco OB2 association, 
even in this case, a scenario s imilar to that described here 
may occur. IWhitworthl lll979T) has shown that the expan- 
sion of an ionised region into molecular gas is preceded by 
a dense shocked layer in the molecular region. In this case 
however, the properties of clumps would probably be differ- 
ent, with those in the shocked layer (and in this case, the 
moving population) begin denser than those in the rest of 
the molecular cloud. Although we have demonstrated in this 
paper that the exact form of the flows has no real influence 
on the resulting clump mass spectra, we have only examined 
the simplest case of initially identical clumps. Further inves- 
tigation would be required to establish whether a triggering 
event from a nearby OB association could produce the effect 
we describe in this paper. 

We note here that our simulations do not contain any 
model for magnetic fields, and as such, we can say nothing 
about such fields wou l d affe ct the clump mass distribution. 
IClifford fc Elmegreenl lll983t) have shown that magnetic link- 
ing between clumps can increase the chance of clump-clump 
interactions. At low Mach numbers, these magnetic interac- 
tions may be stronger than direct clump-clump collisions, 
thus dominating the processes outlined in this paper. 



6 CONCLUSIONS 

We have demonstrated that the process of colliding clumpy 
flows can result in a clump mass distribution with a power 
law consistent with Salpeter (7 = 2.35). Our results high- 
light th at the clump mass distribution observ ed in the 
o Qph (IMotte et alJll998l: I Johnstone et al|l2000|). Serpens 
llTesti ^^argentT^9aT MLd the Orion B Oohnstone et alJ 
l20tnF ~molecular clouds, need have nothing to do with the 
on going star formation. Instead, clumpy structure may arise 
naturally if these clouds were formed from large scale clumpy 
shocks. 

Our results are largely indifferent to our choice of flow 
Mach numbers and filling factors, or even whether gravity is 
included in the model, or not. This suggests that the struc- 
ture and kinematics of the gas which forms star forming re- 
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Figure 4. In this figure we show the clump mass spectra for the 40% filling factor flows but this time, the self gravity of the gas is 
included in the simulations. The un-shaded region denotes the full clump distribution. The hatched region is the sub set of clumps for 
which \E grav \ > 

^thermal and the black regions denote the clumps for which |7£gy olI | ^ (^thermal H~ ^kinetic)- Again, the clump spectra 
are plotted at 0.5£ C r and t cv . The long dashed line shows a Salpeter slope of y = 2.35 (so —1.35 on our log-log plot here). 
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gions, is not important in forming the observed clump/core 
mass distribution. 

The mass spectra in our simulations are not however 
constant in time. Only at the densest point of the shock 
process (i.e. when the full length of the flows have entered 
the shock region) do all the mass spectra yield the Salpeter 
type slope. This is however, the point at which our self- 
gravitating simulations show us that star formation sets in, 
and the spectra exist for the embedded phase: the time pe- 
riod in which we are likely to try and carry out observations 
of clump/core populations. 
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